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ABSTRACT 

The  objective  of  the  ASL  research  effort  in  acoustic  propagation  is 
to  provide  the  Army  with  a  multi -stream  model  for  investigating 
acoustic  detection  systems.  The  first  step  in  developing  this  model 
is  to  account  for  turbulent  scattering.  Five  elements  are  neces¬ 
sary  to  accomplish  this  step--(l)  model  the  turbulent  region  as  a 
collection  of  vortices  with  a  distribution  of  characteristic  sizes/ 
velocities;  (2)  characterize  each  vortex  (turbule)  as  a  known  (or 
assumed)  velocity  distribution  in  three  space;  (3)  solve  the  fluid 
equations  to  determine  the  scattering  from  each  turbule;  (4)  sum  the 
contributions  to  the  scattered  sound  pressure  level  at  the  detector 
location  of  all  turbules  accounting  for  the  propagation  characteris¬ 
tics  of  the  atmospheric  medium;  and  (5)  incorporate  the  algorithms 
devised  above  into  existing  (or  appropriately  modified)  propagation 
models.  Progress  in  these  five  areas  will  be  reported. 

1.  INTRODUCTION 

Acoustic  propagation  in  the  atmosphere  is  influenced  by  a  number  of  factors 
that  have  to  be  considered  in  any  model  designed  to  predict  detector  response 
at  battlefield  distances.  Figure  1  illustrates  a  possible  scenario  that,  along 
with  fig.  2,  incorporates  most  of  the  factors  of  interest.  Figure  1  shows  an 
object  of  interest  marked  S  that  is  a  source  of  the  acoustic  field  marked  F. 
An  acoustic  detector  marked  D  is  some  distance  from  the  source.  An  obstruction 
marked  0  illustrated  as  a  ridge  intervenes  between  S  and  D.  Atmospheric  turbu¬ 
lence  is  marked  T  and  the  reflecting  ground  surface  is  marked  P.  E,  M,  and  R 
are  the  diffracted,  turbulence  scattered,  and  reflected  fields,  respectively. 
The  electronic  signal  from  D  marked  V  is  produced  by  a  superposition  of  the 
direct  field  F  (if  any)  with  E,  M,  and  R  plus  the  intrinsic  detector  noise 
denoted  by  N.  The  problem  of  interest  is  to  determine  theoretically  the  total 
field  given  the  scenario.  The  importance  of  acoustics  as  a  non- line -of -sight 
detection  tool  for  the  Army  is  also  illustrated  in  fig.  1.  The  detector  is 
shown  in  what  is  known  as  a  shadow  zone  by  analogy  with  ray  propagation  from 
optics . 

A  shadow  zone  can  also  be  produced  by  an  upwardly  refracting  atmosphere.  This 
effect  is  illustrated  in  fig.  2.  In  the  figure,  a  source  is  located  at  an  ele¬ 
vation  of  10  m  in  an  atmosphere  with  a  constant  sound  speed  gradient.  The  ray 


‘formerly  U.S.  Army  Atmospheric  Sciences  Laboratory  (ASL) 


43 


path  shown  just  grazes  the  ground  and  is  known  as  the  limit  ray.  The 
shadowzone  is  to  the  right  of  the  curve  minimum  and  between  the  curve  and  the 
range  axis.  The  indicated  height  at  a  range  of  7000  m  is  174.89  m.  The 
shadow  zone  begins  at  a  range  of  1351.32  m.  Current  acoustic  propagation 
models  do  not  account  for  terrain  induced  shadow  zones  as  depicted  in  fig.  1 
but  are  able  to  model  propagation  in  a  stratified  atmosphere  whose 
characteristics  vary  only  with  height  as  that  depicted  in  fig.  2. 


Figure  1.  Acoustic  detection  Figure  2.  Limit  ray  coordinates 

scenario.  locating  the  shadow  zone. 


2.  WAVE  SOLUTIONS  FOR  ACOUSTIC  FIELDS 

While  the  ray  solution  shown  in  fig.  2  accounts  for  the  constant  atmospheric 
sound  speed  gradient,  it  does  not  account  for  the  reflections  that  occur  at  the 
ground  nor  does  it  account  for  the  wave  nature  of  sound.  Figures  3  and  4 
illustrate  the  effect  of  the  ground  but  not  the  full -wave  nature  of  sound.  As 
a  check  on  the  more  sophisticated  models,  a  propagation  model  called  here 
SPHERICAL  was  written  that  models  a  uniform  atmosphere  situated  above  a  flat 
earth,  the  latter  modeled  by  an  impedance  change.  These  plots  using  data 
produced  by  SPHERICAL  are  essentially  ray  solutions  with  the  ground  effect 
modeled  as  an  image  source  below  ground  level,  with  the  image  field  phase  and 
amplitude  modified  according  to  the  impedance  discontinuity.  Figure  3  shows 
a  contour  plot  from  SPHERICAL  that  represents  the  sound  field  for  a  uniform 
atmosphere,  frequency  of  170  hertz,  source  height  of  10.0  m,  sound  speed  335.1 
m/s,  and  a  ground  impedance  ratio  of  8.44  +  J10.0.  The  heavy  line  is  the  limit 
ray  position  from  fig.  2.  Figure  4  shows  the  region  near  the  origin  in  greater 
detail  with  the  horizontal  and  vertical  scales  the  same.  Interference  between 
the  direct  and  the  ground  reflected  waves  is  clearly  demonstrated. 

To  demonstrate  the  full-wave  nature  of  sound  propagation,  a  more  sophisticated 
model  is  necessary.  The  fast  field  program  (FFP)  (Raspet,  1985)  is  a  full-wave 
solution  to  the  propagation  equation  that  accounts  for  a  nonuniform  atmosphere 
using  horizontal  planar  layers,  the  propagation  constants  of  each  layer  being 
uniform.  Ground  effects  are  introduced  by  an  impedance  boundary  just  as  in  the 
SPHERICAL  program.  Since  the  ultimate  goal  was  to  model  turbulence  effects  as 
scattering,  FFP  was  included  as  a  subroutine  in  a  driver  that  produced  data 
from  multiple  sources,  multiple  scatterers,  and  multiple  detectors.  This  pro¬ 
gram  is  named  the  acoustical  multistream  propagation  program  (AMPP)  (Auvermann, 
1992).  AMPP  loops  through  a  set  of  sources,  a  set  of  scatterers,  and  a  set  of 
detector  positions,  saving  the  sound  field  results  of  each  calculation.  The 
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results  for  the  detector  position  set  provides  the  data  for  contour  plots  shown 
below.  AMPP  preserves  the  sound  pressure  amplitude  and  phase,  as  well  as  the 
exitance  (watts/square  meter)  so  that  a  coherent  summation  is  possible  for  a 
reasonable  number  of  sources  or  scatterers.  This  coherent  summation  capability 
of  AMPP  has  not  been  exercised  because  of  the  diversion  caused  by  the 
interesting  feature  of  the  field  in  shadow  zones  that  is  discussed  below. 

To  exercise  the  AMPP  program,  the  scenario  of  fig.  1  was  calculated  for  an 
upward  refracting  atmosphere.  The  results  are  shown  in  figs.  5  and  6.  The 
uniform  sound  speed  gradient  was  -3.67  •  10-3  s_1,  with  speed  at  zero  height  of 
335.1  m/s  and  a  ground  impedance  ratio  as  in  the  SPHERICAL  program  results  of 
figs.  3  and  4.  The  heavy  line  in  each  figure  is  the  limit  ray  path  from  figure 
2.  As  seen  by  comparison  of  fig.  5  with  fig.  3,  SPHERICAL  and  AMPP  show 
virtually  identical  contours  out  to  a  range  of  1000  m.  Beyond  1000  m,  the  AMPP 
upper  contour  parallels  the  limit  ray  path  to  a  high  degree  of  accuracy.  The 
plot  in  fig.  6  shows  in  greater  detail  the  sound  exitance  level  in  the  shadow 
zone  below  the  limit  ray  path  for  the  range  interval  5000  m  to  7000  m. 

Finding  unexpectedly  high  sound  levels  in  the  shadow  zone  has  diverted  the 
investigation  of  turbulence  scattering  into  this  region  because  the  expected 
scattered  level  was  approximately  the  same.  The  nature  of  the  nonturbulence 
related  sound  level  was  investigated.  To  eliminate  the  possibility  that 
computational  precision  was  inadequate,  AMPP  was  changed  to  64-bit  precision. 
The  code  was  run  at  higher  frequencies,  which  should  produce  faster  roll-off 
at  the  shadow  zone  boundary.  Changes  of  ground  impedance  by  an  order  of 
magnitude  in  either  direction  and  a  three-to-one  increase  of  sound  speed 
gradient  were  tried.  All  these  changes  produced  no  material  change  in  the 
results.  The  conclusion  was  that  the  sound  levels  reported  by  AMPP  are  real 
and  to  be  expected.  The  mechanism  is  most  likely  the  reflection  of  a  spherical 
wave  by  a  plane  boundary,  which  will  produce  a  field  somewhat  more  complicated 
than  that  of  the  reflection  of  a  plane  wave  at  a  plane  boundary.  To  expedite 
the  turbulence  investigation,  an  artificially  attenuated  shadow  zone  contour 
picture  was  constructed  as  reported  in  section  5 . 


Figure  3 .  Uniform  atmosphere 
contours . 


SPHCaiCftL  PB06BAM  C0NT0UBS  NCMJ  THE  S0UBCE 


Figure  4.  SPHERICAL  program 

contours  near  the  source . 
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Figure  5.  Upward  refracting 

Figure  6.  Expanded  AMPP  program 

atmosphere  results.  results  for  longer  ranges. 


3.  SCATTERING  BY  TURBULENCE 

Following  the  "structural  approach"  to  turbulence,  acoustical  scattering  from 
turbulence  is  by  analogy  with  optical  scattering  from  aerosols  considered  as 
scattering  from  a  collection  of  eddies  (or  turbules  as  hereinafter  referred  to) 
of  different  sizes  with  a  number  concentration  assigned  to  each  size.  If  the 
scattering  cross  section  and  number  concentration  for  each  size  class  is  avail¬ 
able,  then  a  scattering  calculation  can  be  made  using  the  geometry  of  fig.  7. 
Figure  7  shows  a  source,  a  detector,  and  a  turbule  positioned  over  flat  ground. 
The  coordinate  system  to  be  used  is  also  shown.  The  horizontal  coordinate 
variable  is  r  and  the  vertical  coordinate  variable  is  z.  The  scattering  angle 
variable  6  is  also  shown.  The  quantities  Rat,  Rtd,  and  Rad  are  slant  ranges. 
In  this  paper,  all  distances  are  expressed  in  meters.  The  coordinates  of  the 
principal  objects  chosen  for  this  scenario  and  the  associated  slant  ranges  are 
as  follows : 

ra  -  0.0,  za  -  10.0  R#t  -  5,053.588 

rt  -  5000.0,  zt  -  744.0  Rtd  -  1,246.409 

rd  -  6000.0,  zd  -  0.0  Rad  -  6,000.008 

Subscripts  s,  t,  and  d  stand  for  source,  turbule,  and  detector,  respectively. 
With  this  geometric  arrangement  the  scattering  angle  is  6  -  jt/4  radians. 
Symbols  to  be  used  and  not  already  defined  are  listed  in  table  1  along  with 
units  and  values . 

The  value  of  fld  in  table  1  was  chosen  as  representative,  leading  to  the 
determination  of  the  scattering  volume  from  eq.  (1). 

Qd  -  2w[l  -  cos (*)]  rv  *  Rtd  sin(*)  Vaca  -  *  rv2  1  (1) 

Assuming  an  isotropic  source  in  an  infinite,  lossless  uniform  medium,  the 
exitances  at  the  turbule  and  at  the  detector  are 
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EIt  -  E.  (Rg/R.t)2  -  3. 9156* lO-8  LIt  -  10  log10(EIt/Es)  -  -74.07  dB 
Eid  ”  E„  (R^/R^d)2  -  2.7777.10'8  LId  -  10  log10(EId/E,)  -  -75.56  dB 

The  logarithmic  equivalent  of  these  quantities,  called  level  and  expressed  in 
decibels  (dB) ,  is  given  above  also.  The  letter  I  in  the  subscripts  indicates 
these  quantities  are  infinite  space  values.  In  table  1,  Et  was  set  equal  to 
EIt.  Equation  (2)  is  used  to  calculate  the  scattered  exitance. 


osca(a,v,B)  =  na2(i?(1)  (ie,f)>  oaca(v,Q)  =  fdaoaca(a,v,6)  na(a) 


Ed  - 


_  gtogca(v,6)  va, 


(2) 


Rid 


The  nature  of  the  turbule  concentration  function  na(a)  in  eq.  (2)  is  not  known 
at  this  time  so  an  approximate  evaluation  of  the  integral  over  turbule  size 
will  be  made.  This  approximation  will  involve  selection  of  a  typical  size  at 
and  a  corresponding  typical  number  concentration  with  the  integral  being 
given  by  the  product  of  and  osca(at ,  vt ,  ) ,  where  vfc  is  the  velocity  of  the 
typical  turbule.  The  scattering  volume  VBCa  is  that  of  a  length  1  cylinder 
with  a  radius  at  midpoint  set  by  the  detector  field  of  view  fid. 

The  scattering  efficiency  of  turbules  was  determined  previously  (Goedecke, 
1992).  Theory  was  developed  for  acoustical  scattering  by  localized  quasi¬ 
static  atmospheric  turbules  that  contain  both  flow  and  sound  speed  variations. 
Differential  and  total  cross-section  expressions  have  been  found  for  a  Gaussian 
spinning  turbule  model  of  effective  radius  £  by  using  the  first  Born  approxima¬ 
tion,  which  is  estimated  to  be  valid  for  ka  <  5  for  vmax/c«  <0.1.  The  contri¬ 
butions  of  the  sound  speed  variation  are  about  three  orders  of  magnitude 
smaller  than  those  of  the  flow.  The  latter  are  proportional  to  (ka)6  and 
[sin(0)  co  s(0)  cos  ( 0a)  sin(<£  -  <f>a)]2t  where  (6  $<f>)  are  the  (polar,  azimuthal) 
scattering  angles,  (6a,<f>a)  are  the  spin  axis  angles.  The  contributions  of  the 
sound  speed  inhomogeneities  vary  as  (ka)*  and  have  a  somewhat  Rayleigh -like 
dependence  on  6 . 

Analytical  results  from  this  theory  for  turbule  scattering  efficiency  averaged 
over  turbule  spin  axis  orientation  are  used  above.  These  results  are 
summarized  here.  A  model  spinning  turbule  was  chosen  to  have  the  following 
(quasi)  static  flow  velocity  ?(.?)  : 


?(f)  =  (Qxf)exp (-r2/a2)  Q  =  (2  e)1/2vt/a  .  (3) 

This  equation  represents  a  nonuniformly  spinning  turbule,  with  angular  velocity 
parameter  G  about  an  axis  6  through  the  origin  of  coordinates .  The  quantity 
c«  is  the  background  wave  speed  far  from  the  turbule;  the  parameter  a  is  an 
"effective  radius"  of  the  turbule;  and  f  is  the  position  vector.  The  expres¬ 
sion  for  G  relates  the  angular  velocity  to  the  maximum  turbule  velocity  which 
occurs  at  r  -  a/21/2. 

If  an  acoustic  plane  wave  with  propagation  vector  £,  where  k  -  2n/\t  is  inci¬ 
dent  upon  this  model  turbule,  then  the  first  order  in  G  differential  scattering 
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efficiency,  which  is  defined  to  be  the  scattering  cross  section  divided  by  the 
physical  cross-section  ira2,  averaged  over  all  spin  axis  orientations  is 


(Q{1)  ( k,£ )> 


( ka )*  (sin0cos0)2exp[-(J:a)2  (l-cos0)  ] 


.  (4) 


Figure  7.  Turbulence  scattering  geometry. 
TABLE  1.  SYMBOL  DEFINITIONS 


Symbol 

Definition 

Value 

Units 

a 

turbule  characteristic  size 

m 

at 

typical  turbule  size 

17 . 5/* 

m 

V 

turbule  characteristic  velocity 

m^s~1 

vt 

typical  turbule  maximum  velocity 

1.340 

m* s“* 

n 

turbule  angular  velocity 

0.561 

rad®s-1 

A 

acoustic  wavelength 

10.0 

m 

8 

propagation  to  detector  angle 

x/4 

radian 

(0(1)  (Jc,f)) 

differential  scattering  efficiency 

5 . 52*10-5 

asc  «(a,v,tf) 

differential  scattering  cross 
section 

1.49-10"3 

m2 

Rs 

reference  distance  from  source 

1.0 

m 

E, 

exitance  of  source  at  Rs 

1.0 

w«nT2 

Et 

exitance  at  scattering  volume 

3.92-10"8 

w«nf2 

Ed 

scattered  exitance  at  detector 

9 . 99*10-12 

w»nT2 

«d 

detector  solid  angle  field 

1.0 

str . 

n,(a) 

turbule  concentration  distribution 

Df3 

typical  turbule  concentration 

(1/17. 5)3 

m”  3 

1 

length  of  scattering  volume 

1000.0 

m 

radius  of  scattering  volume 

674.6 

m 

detector  field  half- angle 

0.572 

radian 

V.c. 

scattering  volume 

1.43-109 

m3 
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Figures  8  and  9  show  the  turbule  scattering  properties  to  first  order  in  Oa/c.. 


Figure  8.  Average  differential  Figure  9.  Total  average  scattering 

scattering  efficiency.  for  a  spinning  turbule. 


The  typical  value  for  turbule  size  was  chosen  from  the  size  parameter  for  which 
the  cross  section  is  near  maximum.  Solving  for  at  gives  the  value  in  table  1. 

2*  at/A  -  3.5  <5> 

Since  the  diameter  of  the  turbule  is  approximately  35/w,  a  cell  size  is  chosen 
to  be  17.5  or  slightly  larger  than  the  turbule.  The  cell  is  thought  of  as  a 
cube  that  completely  contains  one  turbule .  The  number  concentration  will  be 
the  reciprocal  of  the  cell  volume,  which  is  the  value  for  nt  given  in  table  1. 
The  final  parameter  to  be  determined  is  the  typical  turbule  velocity.  To  make 
a  reasonable  choice  of  vt,  it  is  appropriate  to  choose  an  outer  scale  size 
related  to  the  dimensions  of  the  scattering  volume  and  to  choose  the  velocity 
of  the  maximum  sized  turbule  as  a  reasonable  wind  speed.  Thus,  if  the  maximum 
turbule  size  is  considered  to  be  500.0  m  and  the  wind  speed  is  6.0  m  •  s  ,  the 
one-third  power  law  of  velocities  determines  the  typical  maximum  velocity  from 
the  relation 

vt  -  6.0  (at/500.0)1/3  ,  <6) 

again  giving  the  value  in  table  1.  The  differential  scattering  efficiency  is 
determined  by  using  eq.  (4);  the  turbule  scattering  cross  section,  the 
scattering  cross  section  per  unit  volume,  and  the  scattered  exitance  at  the 
detector  are  determined  by  using  eq.  (2). 

4.  COMPARISON  OF  RESULTS 

The  purpose  of  this  section  is  to  compare  the  scattered  exitance  derived  in 
section  3  with  that  calculated  from  the  "statistical  approach"  relations  as 
well  as  with  experimental  results  (Gilbert,  1990) .  The  standard  statistical 
approach  scattering  cross  section  is  (Tatarskii,  1971) 
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oTat(0)  =  0.38  k1^ 


(cos  (6) ) 2 


(2 sin(0/2) ) 11/3 


if  (cos (0/2)) 2  +0.13  — 

cl  it 


(7) 


The  symbols  that  have  not  been  defined  before  are  Cv2,  the  velocity  structure 
constant,  CT2,  the  temperature  structure  constant,  and  T®,  the  ambient 
temperature.  Equation  (7)  has  been  rederived  with  fewer  approximations  but  the 
result  has  not  been  reported.  The  derivation  is  too  lengthy  to  include  here, 
but  the  final  result  is 


=  2tiJc4 


ElUsl.  (cos  (0)  cos  (0/2)) 2  + 

4i i 


cl 


Bv(k)  =4.16  C*  BT(k)  =  4 . 16  Cy 

cl  =  0.04  (in  •  s-x)2m-2/*  c£  =  0.00043  ic2in '2/3  . 


(8) 


The  numerical  constant  4.16  (dimensions  of  m11/3)  in  the  expressions  for  Bv(k) 
and  BT(k)  in  eq.  (8)  follows  approximately  the  angle  dependence  of  the  first 
bracket  in  eq.  (7)  and  has  been  determined  for  an  angle  of  x/4  radians.  The 
structure  constants  were  taken  from  the  literature  of  Brown  (1976).  When  an 
asymptotic  wave  speed  of  335.1,  a  temperature  of  300  k,  and  eqs .  (7)  and  (8) 
are  used,  the  statistical  and  Tatarskii  cross  sections  are 

a8ta(0)  -  6.2417  •  1CT7  m"1  aTat(0)  -  1.3217  •  1CT7  nf1  .  (9) 

Detected  exitance  using  eq.  (2)  and  excess  loss  levels  for  the  four  cases  of 
interest  are  brought  together  in  table  2 . 

TABLE  2.  DETECTED  FIELDS  COMPARISON 


Symbol 

Definition 

Value 

Units 

Ed 

exitance,  structural  approach 

9 . 99*10“ 12 

w«m~2 

E.ta 

exitance,  statistical  approach 

2 . 25*10-11 

w*nT2 

^Tat 

exitance,  Tatarskii  formula 

4. 76«10-12 

w»m-2 

Ld 

excess  loss,  structural  approach 

34.44 

dB 

f*sta 

excess  loss,  statistical  approach 

30.91 

dB 

f^at 

excess  loss,  Tatarskii  formula 

37.66 

dB 

Llit 

excess  loss,  literature  value 

20  -  30 

dB 

A  number  of  issues  impact  the  validity  of  the  excess  loss  values  in  table  2. 
L^  was  calculated  using  an  assumed  turbule  number  concentration  and  velocity  and 
a  typical  turbule  size  rather  than  an  integration  over  a  size  distribution. 
Lata  was  calculated  using  eq.  (8),  which  was  derived  using  a  realistic  structure 
function  that  became  constant  at  large  seperations .  Perhaps  this  is  the  reason 
for  the  difference  between  Lata  and  Lrat.  and  Lata  are  near  the  spread  in  Lut. 
The  latter,  however,  was  measured  for  a  different  geometry  and  different 
frequency.  The  fact  that  the  new  results,  as  crude  as  represented  here,  give 
answers  in  general  agreement  with  the  statistical  approach  and  with  the 
experiment  is  encouraging. 
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5.  TURBULENT  SCATTERING  IN  A  PROPAGATION  MODEL 


Figures  10  through  12  show  inclusion  of  turbulent  scattering  in  a  propagation 
model  by  running  AMPP  twice  and  summing  the  results.  The  first  run  preserves 
the  field  produced  by  the  source  (fig.  10).  Figure  10  is  a  reproduction  of 
figure  6  except  the  shadow  zone  has  been  artificially  cleared  of  residual  sound 
levels.  The  second  run  preserves  the  field  produced  by  an  isotropic  scatterer, 
the  results  being  shown  in  fig.  11.  The  two  results  are  combined  after  the 
second  run  results  have  been  modified  by  the  source  to  turbule  loss.  This 
final  result  is  shown  in  fig.  12  where  the  summation  has  been  done 
incoherently.  Figure  12  only  gives  a  representative  picture  of  the 
superimposed  fields  because  the  solution  algorithm  of  AMPP  (from  FFP)  can  only 
model  isotropic  scatterers;  whereas,  it  is  known  that  scatterers  have  patterns 
resembling  fig.  8.  The  scattered  exitance  in  fig.  12  has  been  adjusted  to  give 
an  excess  loss  at  the  detector  location  of  fig.  7  of  approximately  25  dB. 


6. 


CONCLUSIONS 


The  "structural  approach"  to  modeling  turbulence  produces  scattered  sound 
levels  that  are  in  general  accordance  with  those  calculated  using  the 
"statistical  approach."  They  are  also  in  accordance  with  the  experiment.  The 
claim  is  not  that  a  turbulent  region  is  made  up  only  of  a  collection  of 
turbules  of  various  sizes.  Rather,  the  claim  is  that  such  a  concept  is  a  first 
approximation  to  the  true  distribution  of  motion,  energy,  and  density  in  the 
turbulent  region.  Furthermore,  as  far  as  scattering  is  concerned,  this  concept 
yields  reasonable  answers.  The  significance  of  the  work  reported  here  is  that 
a  pathway  has  been  charted  that  breaks  the  shackles  imposed  by  the  isotropic 
and  homogeneous  requirements  of  the  statistical  approach.  Boundary  layer 
theory  and  the  vertical  structure  capability  of  current  acoustic  propagation 
models  both  attest  to  the  recognized  fact  that  the  atmosphere  is  not  isotropic 
and  homogeneous.  The  structural  approach  supported  by  analytical  or  numerical 
knowledge  of  the  scattering  properties  of  turbules  does  work.  What  has  not 
been  done  yet  is  to  determine  the  number  concentration  distribution  appropriate 
for  use.  Additionally,  knowledge  of  the  scattering  properties  of  large 
turbules  is  yet  to  be  determined.  The  true  internal  structure  of  a  turbule  has 
not  been  determined,  but  the  scattering  results  do  not  seem  to  be  very 
sensitive  to  this  structure.  Use  of  the  structural  approach  requires  that  the 
analyst  know  something  more  of  the  details  of  the  turbulent  region  than  the 
statistical  approach.  A  simulation  of  flow  in  the  actual  environment  being 
analyzed  is  indicated.  Calculation  of  scattering  from  such  unusual  features 
as  intermittent  turbulent  regions  is  possible.  It  will  be  necessary  to  devise 
a  method  by  which  anisotropic  sources  and  scatters  can  be  included  in 
propagation  models . 
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